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ABSTRACT 

Objective: Adherence with bisphosphonates therapy 
is generally low. Enhancing adherence with 
bisphosphonates would be effective in achieving the 
full benefits of therapy albeit a growth in the 
expenditure for supporting incremented drug use is 
expected. The cost-effectiveness of enhancing 
adherence with oral bisphosphonates in a large 
population of osteoporotic women has been assessed 
in the current study. 
Design: Retrospective cohort study. 
Setting: Healthcare utilisation databases of Lombardy 
Region, Italy. 

Participants: A cohort of 28 558 women aged 
45 years or more, resident in the Italian Region of 
Lombardy, who were newly treated with oral 
bisphosphonates during 2003-2004, was followed for 
6 years after index prescription. 
Outcome measures: Fracture-free survival time, 
healthcare cost and incremental cost-effectiveness ratio 
(ICER) of enhancing adherence, that is, the additional 
cost that would be spent every year for gaining one 
fracture-free year as a consequence of enhancing 
adherence at a certain level. 
Results: Enhanced adherence from 33% (baseline) to 
80%, increased both fracture-free survivals from 970 
to 973 years and healthcare costs from €118 000 to 
€265 000 every 1000 woman-years, with ICER value of 
€53 000 (95% CI €49 000 to €58 000). ICER values 
were lower for older women (€50 000; 95% CI 
€42 000 to €58 000) and for those suffering from at 
least a chronic comorbidity (€25000; 95% CI 95% CI 
€7000 to €47 000). 

Conclusions: Enhancing adherence with oral 
bisphosphonates offers important benefits in reducing 
the risk of fracture, although at a substantial cost. 
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INTRODUCTION 

Osteoporosis is a silent systemic disease char- 
acterised by low bone mineral density 
(BMD) and alterations of the skeleton's 
microarchitecture that determine fragility of 



Strengths and limitations of this study 



The study was conducted in a large unselected 
population reflecting the patients' treatment in 
the real-life setting. 

A peculiarity of the current approach is that esti- 
mates are entirely based on observed data in the 
setting of daily clinical practice, rather than on 
usual pharmacoeconomic analysis. 
Residual confounding might partially affect the 
study results. 



the bone and subsequent increased risk of 
fracture even in the case of mild traumas. 1 
With ageing populations, the burden of 
osteoporotic fractures on society will increase 
in the coming years 2 and the prevention of 
osteoporotic fractures is therefore a major 
public health issue. 

Randomised clinical trials (RCTs) have con- 
sistently shown that long-term treatment with 
oral bisphosphonates (BPs) is effective in 
improving BMD from 2.4% to 8.8% and redu- 
cing the risk of fracture from 12% to 49%. 3 4 
However, the efficacy found in RCTs may not 
be replicated in actual clinical practice, 
mainly because patients enrolled in RCTs 
achieve almost optimal treatment adherence, 
while in the clinical practice adherence is 
impaired. Poor adherence with BPs therapy 
has been observed in the real-life setting. 5 6 
Evidence from observational studies shows 
that poor adherence is associated with sub- 
optimal improvement of bone density, 7 
increasing risk of fracture 5-8 and worsening 
cost-effectiveness profile of therapy. 9 

Enhanced adherence with BPs is thereby 
expected to be effective in achieving full ben- 
efits from therapy. 10 However, increasing 
expenditures for supporting increased drug 
use are also expected as adherence improves. 
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The expenditures' impact is particularly important from 
a public health perspective, since it is important to opti- 
mally allocate the limited available economic resources 
as to maximise the population health level. 11 It would 
therefore be suitable to evaluate cost and effectiveness 
of treatment jointly, in order to quantify the additional 
cost needed to increase the treatment's effectiveness by 
enhancing adherence. 

Model-based economic evaluations are commonly used 
to compare alternative treatment strategies in osteopor- 
osis, to support decision makers and to inform treatment 
guidelines. 12-14 A large part of these approaches is based 
on data on treatment-efficacy obtained from RCTs. 
Although economic analyses carried out alongside RCTs 
have many methodological advantages, their findings 
may be poorly generalisable to the clinical practice. 15 
Other than achieving almost optimal compliance, 5 6 
RCTs often impose strict inclusion and exclusion criteria, 
while patients presenting to physicians for treatment fre- 
quently exhibit features that would exclude them from 
entering a trial. 16 Therefore, observational studies based 
on high-quality data reflecting use and impact of drugs in 
the real world clinical practice might be alternative data 
sources useful for economic healthcare evaluations. 17 

We developed an approach to estimate the cost- 
effectiveness profile of enhancing adherence with BPs 
therapy in a large population without signs of previous 
fracture. Costs and effectiveness were estimated from a 
population-based cohort study. Data were derived from 
the healthcare utilisation (HCU) databases of the 
Lombardy Region (Italy). 

METHODS 
Setting 

The data used in the present study were retrieved from 
the HCU databases of Lombardy, a Region of Italy 
which accounts for about 16% (9 million) of its popula- 
tion. In Italy, the population is fully covered by the 
National Health Service (NHS), which in Lombardy has 
been associated since 1997 with an automated system of 
databases collecting a variety of information, including: 
(1) an archive of residents receiving NHS assistance 
(practically the whole resident population), reporting 
demographic and administrative data; (2) a database on 
discharge diagnoses from public or private hospitals and 
(3) a database on outpatient drug prescriptions reim- 
bursable by the NHS. 

For each patient we linked the above databases via a 
single identification code. In order to preserve privacy, 
each identification code was automatically converted to 
an anonymous code. The inverse process was prevented 
by the deletion of the conversion table. Full details of 
the procedure have been reported elsewhere. 18 

Cohort selection and follow-up 

The cohort regarded all women aged 45 years or older 
who received at least one prescription of oral BPs 



(alendronate or risedronate) between 1 January 2003 and 
31 December 2004. According to current Italian regula- 
tions, these women certainly received diagnosis of osteo- 
porosis, since BPs therapy reimbursement is obtained only 
for prophylaxis of osteoporotic fractures. The first pre- 
scription of BPs was defined as the index prescription. 

To make the data as relevant as possible to the study 
aim, five categories of patients were excluded. One, 
patients who within 3 years before the index prescription 
had received BPs drug therapy, in order to ensure the 
inclusion of only newly treated women. Two, patients who 
within 3 years before the index prescription, had received 
drugs known to reduce the risk of bone fracture (ie, calci- 
tonine, raloxifene and/ or hormone replacement therapy 
(HRT)). Three, patients who within 3 years before the 
index prescription, had been hospitalised for the 
outcome of interest, to focus the data on primary preven- 
tion of fracture. Four, patients who within 3 years before 
the index prescription, had been hospitalised for diseases 
known to increase the risk of bone fracture (ie, cancer, 
Paget's disease, diabetes, thyroid disease, chronic renal 
failure and/or disease of connective tissue). Five, patients 
who did not reach at least 1 year of follow-up, to ensure at 
least 1 year of potential exposure at the drugs of interest. 

Each cohort member accumulated person-years of 
follow-up from the date of index prescription until the 
earliest among the dates of outcome occurrence 
(ie, hospital admission for fracture) or censoring, that is, 
death, emigration or end of follow-up (fixed at 6 years 
after the date of index prescription). 

Exposure assessment 

Drug type, dosage, regimen and dispensation date of all 
BPs prescribed to each cohort member during follow-up 
were retrieved from the outpatient drug prescriptions 
database. During the study period two drug types (alen- 
dronate and risedronate), either on once daily (10 and 
5mg/day, respectively) or once weekly (70 and 35 mg/ 
week, respectively) regimens, were available for reim- 
bursement by NHS. The time span with drug available 
for each prescription was calculated assuming the stand- 
ard frequency of intake and the prescribed dosing 
regimen (ie, 2 or 4 weeks), and the date of exposure 
ending of each dispensation was established accordingly. 
In this way, each day of the follow-up period was evalu- 
ated as covered or not covered by BPs availability, irre- 
spective of whether the availability regarded single or 
overlapping prescriptions. Adherence was measured by 
dividing the cumulative number of days with BPs avail- 
able for the number of days of follow-up, a measure 
denoted as proportion of days covered (PDC). 19 
A patient was considered adherent to BPs therapy if her 
PDC value was 80% or higher, a value conventionally 
used in most database studies on this topic. 8 

Outcome ascertainment 

Hospitalisations for traumatic or pathological bone frac- 
ture (International Classification of Diseases (ICD)-9 CM: 
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800-829, 733.1) occurred to cohort members during 
follow-up were retrieved from the hospital discharge data- 
base, and the earliest date of hospitalisation was consid- 
ered as that of outcome occurrence. 

Additional measurements 

Several factors measurable in the HCU database of 
Lombardy were considered. These included age at the 
date of index prescription, BP type (alendronate or rise- 
dronate) and regimen (daily or weekly) filled at therapy 
initiation, and other drugs dispensed during follow-up 
(HRT, antidepressants, antidiabetic agents, lipid-lowering 
drugs, diuretics and (3-blockers). In addition, the 
Charlson comorbidity index, 20 was calculated using diag- 
nostic information from inpatient encounters in the 
3 years prior and 1 year after the date of index prescrip- 
tion and expressed in two categories (ie, 0 or 1, respect- 
ively, suggesting the absence or presence of at least one 
comorbidity sign) . 

Estimating the relationship between adherence 
and outcome 

The Cox proportional hazard regression was used to esti- 
mate the HR, and corresponding 95% CI, for the associ- 
ation between adherence with BPs and time of outcome 
occurrence. Estimates were adjusted for time-fixed 
factors (age, BP type and regimen dispensed at therapy 
initiation and Charlson comorbidity index) as well as for 
factors measured during follow-up (adherence with BPs 
therapy and cotreatments with other drugs). As adher- 
ence, as well as all other factors measured during 
follow-up, can change over time, assessment of their 
effects requires properly accounting for the varying 
nature of these variables. This was carried out by fitting 
the Cox model that includes these factors as time- 
dependent covariates. 21 For instance, by considering the 
predictor variable of interest (ie, PDC categories), with 
this approach each participant's cumulative adherence is 
recalculated from the start of follow-up until the time of 
each outcome onset. Thus, the HR associated with a cat- 
egory of PDC is derived using information concurrent to 
the observed outcomes, rather than the adherence 
profile over the full length of follow-up. The same 
approach was used for all other time-varying factors 
above reported. 

Estimating cost and effectiveness of enhancing adherence 

Starting from the observed average PDC (baseline), 
three scenarios were built by progressively enhancing 
the adherence average level at 40%, 60% and 80%. 
Healthcare costs and fracture-free survivals observed at 
baseline, as well as those expected for each given scen- 
ario, were computed as follows. 

First, we considered the observed and expected distri- 
butions of the included patients according to their PDC 
values. For each scenario, a new PDC value was assigned 
to each patient by adding to her baseline PDC a random 
increment drawn from a truncated normal distribution 22 



with mean and variance derived from the difference 
between the average PDC observed at baseline and that 
imposed from the scenario. 

Second, health outcomes were estimated by means of 
the years spent free from fracture observed at baseline, 
as well as those expected for a given scenario. Observed 
fracture-free survival time was computed following 
Neymark et al. 23 Briefly, the Breslow estimator was used to 
calculate the baseline hazard of each cohort member in 
presence of time-dependent covariates. 24 The baseline 
hazard was combined with the parameters estimated 
from the Cox model to derive the survival curve of a 
generic cohort member, characterised by a specific 
history of exposure and covariate values. The years spent 
free from fracture by each cohort member were esti- 
mated as the area under the survival curve 24 up to the 
maximum observed follow-up time. 25 Individual survival 
years were then summed so resulting in the observed 
number of fracture-free years for every 1000 woman- 
years. Expected survival years were computed using the 
same procedure, but substituting the observed PDC 
value with that expected in the given scenario, keeping 
the values of all other covariates unchanged. 

Third, healthcare annual costs observed at baseline, as 
well as that expected for a given scenario, were esti- 
mated. Costs from NHS perspective were calculated by 
the amount that the Regional Health Authority reim- 
bursed to health providers. Costs for supporting BPs 
therapy, as well as those for supporting healthcare in the 
1-year time horizon after index hospitalisation actually 
experienced by cohort members, were regarded as of 
interest for calculating observed cost. With the aim of 
providing an order of costs magnitude, the daily price 
for BPs therapy ranged from €0.35 to €0.65, while costs 
for care of a patient experiencing fracture in the first 
year after index hospitalisation (including direct costs 
for inpatients and outpatients care and assistance) 
ranged from €2147 for wrist fracture to €14 221 for hip 
fracture. The same two sources were regarded for calcu- 
lating expected costs. In particular, those for supporting 
BPs therapy at a given level of enhanced adherence (of 
course expected to be increased with respect to base- 
line) were estimated from the distribution of patients 
over the categories of PDC, assuming that prescription 
profile does not change by increasing adherence. Costs 
for supporting patients' healthcare experiencing frac- 
ture for a given scenario (of course expected to decrease 
by enhancing adherence) were estimated from the 
fracture-free years expected for that scenario. 

Finally, we measured the incremental cost-effectiveness 
ratio (ICER), determined by dividing the differences in 
cost and health outcomes between two strategies, that is, 
observed at baseline versus expected according to a 
given scenario of enhanced adherence. 14 The ICER is 
the additional cost that would be spent every year for 
gaining one fracture-free year as a consequence of 
enhancing adherence at a certain level. The ICER's CIs 
were calculated using a non-parametric bootstrap 
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method based on 1000 resamples/ b The ICER was esti- 
mated for the whole cohort, as well as for strata of age at 
the index prescription and Charlson comorbidity score. 

Sensitivity analyses 

To verify the robustness of our estimates, three sensitivity 
analyses were performed and the corresponding ICER 
values were recalculated. One, alternative PDC thresh- 
olds to those used in the main analysis were considered. 
Two, alternative drug costs to those observed were used 
by substituting them with the cost of the least costly drug 
with the same active principle as that dispensed, as well 
as with the least costly drug among those available in 
Italian market. Three, the impact of unmeasured risk 
factors on adherence-outcome association was evaluated 
using Monte Carlo sensitivity analysis. 27 To account for 
confounders potentially oversizing the adherence- 
outcome association, we considered factors known to 
increase the fracture risk that also negatively affect medi- 
cation adherence. Depression was identified as a pos- 
sible source of such bias, since evidence exists that (1) it 
is a risk factor for the occurrence of bone fracture 28 29 
and (2) it is associated to less frequent use of health ser- 
vices, including poor treatment adherence. 30 According 
to a recent meta-analysis, we assumed that depressed 
patients have a threefold greater risk to be non-adherent 
to BPs therapy compared with patients who did not 
suffer from depression. 30 Moreover, we assumed the frac- 
ture risk to be 1.5 or 2.0 times greater in patients with 
depression compared to those without depression, in 
accordance with data from a recent meta-analysis of pro- 
spective studies examining the risk of osteoporotic frac- 
tures associated with depression, 28 and findings from the 
US National Health and Nutrition Examination 
Survey. 29 Following the procedure described by 
Steenland and Greenland, 2 the adherence-fracture 
association was calculated after adjusting for such 
unmeasured confounder used to recalculate the ICER. 

The SAS statistical package was used for the analyses 
(SAS, V.9.2). All the tested hypotheses were two-tailed, 
with p value <0.05 considered statistically significant. 



RESULTS 

Cohort characteristics 

In total, 57 175 women aged 45 years or older with at 
least a BPs prescription during 2003-2004 were identi- 
fied. Among these, 24 274women were excluded 
because they had already received BPs or other cotreat- 
ments. Further 3273 and 1070 women were, respectively, 
excluded because they were affected by selected dis- 
eases, or had less than 1 year of follow-up. The remain- 
ing 28 558 women were included in the cohort. They 
accumulated 151 382 woman-years of follow-up (on 
average 5.3 years per patient) and 50 428 woman-years 
of therapy (on average 33% of the time of follow-up was 
spent with BPs available). The cohort members who 



during follow-up experienced at least one hospitalisation 
for fracture were 1921. 

Table 1 summarises selected characteristics of cohort 
members according to their adherence with BPs therapy 
during follow-up. Mean age at cohort entry was 72 years 
(SD 9.6 years). Alendronate administered once weekly 
was the most employed first-choice therapy. Overall, only 
14.6% of patients were adherent with BPs therapy (ie, 
>80% of follow-up was covered by drug availability). 
Compared with women who adhered with BPs therapy, 
those with poor adherence were slightly younger (almost 
1 year) , were initially treated with alendronate adminis- 
tered once daily, and were less frequently cotreated with 
corticosteroids and more frequently with antidepressant 
and antidiabetic agents. There was no evidence that 
cotreatment with lipid-lowering drugs, diuretics and 
(3-blockers, as well as that worse profile of comorbidity, 
affected adherence with BPs therapy. 

Relationship between adherence and outcome 

Table 2 shows that women who adhered with BPs 
therapy had a significant 18% reduction in fracture risk 
compared with those who did not adhere to the therapy. 
There was no evidence that either drug type or drug 
regimen administered at index prescription predicted 
the following risk of fracture. The fracture risk was sig- 
nificantly increased among older women who had at 
least a comorbidity. Cotreatment with corticosteroids, 
antidiabetics, lipid-lowering and antihypertensive agents 
were associated with a reduction in fracture risk. 

Cost and effectiveness of enhancing adherence 

Figure 1 shows the trends in annual cost and fracture-free 
survival years, as well as in ICER point and interval esti- 
mates, from baseline to scenarios of progressively 
enhanced adherence with BPs therapy. Progressive 
increase in costs (from €118 000 to €265 000 every 1000 
woman-years) and in the time span free from fracture 
(from 970 to 973 years every 1000 woman-years) was 
observed. With respect to baseline, the annual cost that 
should be additionally spent to gain one fracture-free year 
was €127 000 (95% CI €110 000 to €146 000) when enhan- 
cing adherence at 40%, but fell to €53 000 (95% CI 
€49 000 to €58 000) when enhancing adherence at 80%. 

As expected, ICER values were notably lower among 
older women suffering of comorbidities. The ICERs cal- 
culated for the 80% adherence scenario (compared to 
baseline) were, respectively, €274 000 (95% CI €240 000 
to €310 000), €55 000 (95% CI €36 000 to €77 000) and 
€50 000 (95% CI €42 000 to €58 000) for women aged 
<60 years, from 60 to 70 years and >7l years. The ICER 
values were €72 000 (95% CI €64 000 to €80 000) and 
€26 000 (95% CI €7000 to €47 000) for women without 
and with signs of chronic comorbidities, respectively. 

Sensitivity analyses 

Figure 2 shows that the use of alternative thresholds for 
dichotomising adherence changed to only a small extent 
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Table 1 Selected traits of the included women according to their adherence with bisphosphonates therapy 

Adherence category 



Low coverage* (N 24 383) 



High coverage* (N 4175) 



All patients (N 28 558) 



At the time of index prescription 
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Once weekly 


22 265 (91 .3%) 


3917 (93.8%) 


26 182 (91.7%) 


Cotreatments 








Corticosteroids 


9017 (37.0%) 


1674 (40.1%) 


10 691 (37.4%) 


Antidepressants 


9261 (38.0%) 


1417 (34.0%) 


10 678 (37.4%) 


Drugs used for diabetes 


1777 (7.3%) 


262 (6.3%) 


2039 (7.1%) 


Lipid-lowering agents 


7106 (29.1%) 


1192 (28.6%) 


8299 (29.1%) 


Diuretics/p-blockers 


13 812 (56.7%) 


2398 (57.4%) 


16 210 (56.8%) 


Charlson comorbidity score 








0 


20 176 (82.8%) 


3458 (82.8%) 


23 634 (82.8%) 


>1 


4207 (17.3%) 


717 (17.2%) 


4924 (17.2%) 


Healthcare utilisation database of Lombardy, Italy, 2003-2010. 



*Patients who did not adhere with bisphosphonates therapy were those with a proportion of days of follow-up covered by drug availability 
<80% (low coverage); patients who adhered with bisphosphonates therapy were those with a proportion of days of follow-up covered by drug 
availability >80% (high coverage). 



the cost-effectiveness estimates associated with enhancing 
adherence from baseline to 80% (being the ICER values 
€52 000 and €68 000 for thresholds of 70% and 90%, 
respectively). The ICER values associated with an 80% 
adherence scenario would be reduced to €46 000 or 
€34 000 if all patients were treated with the less costly 
drug within each active principle, or the less costly drug 
among those available in Italian market. Finally, the 
ICERs changed only to a small extent by external adjust- 
ment for depression status, even imposing the strongest 
confounder-outcome association (the ICERs being 
€72 000 and €95 000 when imposing that fracture risk was 



1.5 and 2.0 times greater in exposed than in unexposed 
to the unmeasured confounder, respectively) . 

DISCUSSION 

In our population, nearly 15 000 women aged 45 years 
or older started BPs therapy every year (almost 66 new 
users every 10 000 woman-years). According to both 
Italian regulations and the adopted exclusion criteria, 
these women must have received diagnosis of osteopor- 
osis and must have started therapy for primary preven- 
tion of bone fracture. On an average, the cost for drug 



Table 2 Adjusted HRs, and corresponding 95% CIs, for the association between selected traits of the included women, 



including their adherence with bisphosphonates therapy during follow-up, and the risk of fracture 




Number of fractures 


Adjusted HR* 


95% CI 


Age in years 




1.07 


1.07 to 1.08 


Therapy with bisphosphonates 








At index prescription 








Risedronate vs alendronate 


355 vs 1566 


0.89 


0.78 to 1.00 


Once weekly vs once daily 


1755 vs 166 


0.97 


0.82 to 1.15 


Adherence during follow-up 








High vs low coverage 


344 vs 1577 


0.82 


0.73 to 0.93 


Cotreatments (yes vs no) 








Corticosteroids (users vs no users) 


514 vs 1407 


0.65 


0.58 to 0.72 


Antidepressants (users vs no users) 


744 vs 1177 


1.02 


0.93 to 1.12 


Drugs used for diabetes (users vs no users) 


107 vs 1814 


0.82 


0.67 to 0.99 


Lipid-lowering agents (users vs no users) 


353 vs 1568 


0.59 


0.52 to 0.66 


Diuretics or p-blockers (users vs no users) 


896 vs 1025 


0.72 


0.65 to 0.78 


Comorbidity 








Charlson score (>1 vs 0) 


482 vs 1439 


1.83 


1.65 to 2.04 



Healthcare utilisation database of Lombardy, Italy, 2003-2010. 



*HR estimated by fitting Cox proportional hazard model. Estimates are adjusted for all the factors reported in the table (see text). 
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Figure 1 Trends in healthcare 
cost and fracture-free survival 
(upper box) and in incremental 
cost-effectiveness ratio, with 
corresponding 95% CIs (below 
box), from baseline to scenarios 
of progressively enhancing 
adherence with bisphosphonates 
therapy. Healthcare utilisation 
database of Lombardy, Italy, 
2003-2010. 
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acquisition amounted to about €118 000 every 1000 
woman-years (ie, an annual cost of €1.7 million for the 
entire cohort of BP-therapy initiators). Drug cost, 
however, strongly depends on the average BP-therapy 
adherence level observed in the considered population. 
With respect to other reports, 5 6 a very low adherence 
level was found in our setting, since only 33% of the 
observational period was covered by drug availability. 
Consistently with others, 5-8 we found that, compared 
with women with adherence level <80%, those with 
higher adherence had an 18% reduced risk of fracture. 
This implies that should higher adherence have 
occurred in the population, a lower number of women 
would have experienced hospitalisation for bone frac- 
ture and a gain in fracture free survival time would have 
been observed. Based on our estimates, 0.18 or 2.77 
fracture-free years are, respectively, expected to be 
gained every 1000 woman-years at risk for scenarios 
involving average adherence of 40% or 80% in the 
target population. However, €22 000 (or €146 000) every 
1000 woman-years are expected to be additionally spent 
with these scenarios. Thus, as a balance, we estimated 



that by enhancing average adherence level from baseline 
to 40% (or 80%), €127 000 (or €53 000) would be add- 
itionally spent to gain one fracture-free survival year. 
Additional costs of €51 000 and €26 000 are, respectively, 
expected by enhancing adherence from baseline to 80% 
among older women and among those suffering of 
comorbidities. 

Such estimates should be compared with those sup- 
ported by the NHS for taking care of patients who 
experience bone fracture during the first year after the 
first hospitalisation. Rough estimates may be made by 
dividing the annual direct costs sustained for hospitalisa- 
tions and rehabilitation following bone fractures by the 
number of patients who every year were hospitalised for 
fracture. Estimates of this type may be drawn out by 
recent papers from Italy, France and Canada that, 
respectively, reported annual cost per patient of €11 000, 
€12 000 and €40 000. 31 " 33 

This study is unique in several respects. Results from 
meta-analysis of RCTs have been extensively used to 
model cost-effectiveness of the treatment of osteopor- 
osis. 34 There are, however, circumstances where 
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Figure 2 Influence of choosing PDC (proportion of days covered by bisphosphonates (BPs) availability) threshold for 
dichotomising adherence category (A; see ADHERENCE footnote), spending for drug therapy with BPs (B; see COST footnote), 
and controlling for unmeasured confounders (C; see CONFOUNDING footnote) on the incremental cost-effectiveness ratio 
values associated with enhanced adherence from baseline to 80%. Healthcare utilisation database of Lombardy, Italy, 2003- 
2010. ADHERENCE — in the main analysis a patient was considered adherent with therapy if her PDC value was 80% or higher. 
We reported alternative thresholds of 70% and 90%, respectively, denoted first and second scenario. COST — in the main 
analysis the cost of BPs truly dispensed in the clinical practice was considered. We here considered alternative cost including the 
less costly drug within each therapeutic class as that dispensed (first scenario), and the less costly drug among those available 
in Italian market (second scenario). CONFOUNDING — in the main analysis the extent of adherence-outcome association was 
adjusted for factors directly measured from databases. We here used external adjustment imposing that (1) compared with 
patients unexposed to the unmeasured confounder, those exposed have a three times greater risk to be no adherent with BPs 
therapy; (2) the risk of fracture was 1 .5 (first scenario) and 2.0 (second scenario) times greater in exposed than in unexposed to 
the unmeasured confounder. 



randomisation is not possible, such as in studies aimed 
at investigating the impact of treatment adherence on 
clinical outcomes in the 'real-life' setting. A first peculi- 
arity of our approach is that estimates are entirely based 
on observed data in the setting of daily clinical practice, 
rather than on Markov microsimulation models typical 
of economic analysis. 35 Assuming that a simulation 
model is only necessary when direct observation of the 
phenomenon is not available or cannot be collected 17 
our approach provides an original framework for evalu- 
ating the cost-effectiveness of enhancing adherence 
entirely based on HCU databases. A second peculiarity 
of our approach is that the increase of fracture-free sur- 
vival time was used to measure effectiveness of enhan- 
cing adherence. This approach, although unusual, was 
described as essential for economic evaluations of 
healthcare interventions for which patient survival is the 
principal outcome. 22 Notwithstanding these 



peculiarities, our estimates were of the same order of 
magnitude with respect to those using more usual 
approaches. For instance, a recent microsimulation 
study reported that the costs per QALYs gained, for 
branded BPs were estimated at €19 000, €32 000 and 
€64 000 for medical possession rates of 100%, 80% and 
60%, respectively. 36 

In Italy: (1) a cost-free uniformly organised healthcare 
system involves practically all citizens, (2) the healthcare 
system is managed by a set of databases including inpati- 
ents hospital discharge and outpatients drug prescrip- 
tion data, (3) these databases guarantee universal 
coverage of population receiving health assistance from 
NHS and (4) quality of these databases has been repeat- 
edly ascertained, although in a clinical setting different 
from that considered in the current study. 37 Owing to 
these characteristics, the complete history of health ser- 
vices utilisation of a very large unselected population is 
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available and real world clinical practice, including 
health and economic burden of medical interventions, 
may be investigated. According to our best knowledge, 
the current study offers unique insight on the potentials 
of HCU databases to estimate the cost-effectiveness 
profile of medical interventions on osteoporosis. 

A number of sensitivity analyses allowed us to verify the 
robustness of our findings. For example, the use of alterna- 
tive PDC categorisations as well as adjustment for unmeas- 
ured confounders did not substantially modify the reported 
estimates. We also showed that the cost-effectiveness profile 
was greatly improved by considering less costly drugs, typic- 
ally generic drugs. Of course, this is true under the assump- 
tion that treatment with generic drugs does not lead to 
lower adherence or more fracture-related hospitalisations 
compared with brand-name therapy. 

Our study has a number of potential limitations. 
External validity is a first major concern. Drug costs are 
extremely variable among countries, and will likely 
decrease with increasing competition from generics. 38 
As a consequence, our estimates cannot be universally 
generalised. 

Both the components of cost-effectiveness estimate are 
likely affected by sources of systematic uncertainty. As far 
as costs are concerned, our analysis did not include 
those to support interventions aimed at improving 
adherence. Although interventions that try to change 
the adherence level in population can be very costly and 
the objective is very hard to obtain, 39-41 we were not 
able to measure or even approximate costs for support- 
ing interventions aimed at improving adherence since, 
to our best knowledge, no studies have been performed 
on this topic in Italy. Furthermore, interventions able to 
enhance adherence unlikely have an effect only on one 
therapeutic class, so that their cost should impact on a 
variety of chronic therapies. This implies that our study 
must be regarded as aimed to investigate the cost- 
effectiveness of enhanced adherence, rather than the 
cost-effectiveness profile of interventions aimed at 
improving drug adherence. In the past few years a sub- 
stantial improvement of adherence with chronic drug 
therapies has been noticed in Lombardy Region of Italy, 
though no specific interventions were implemented in 
this regard. This suggests that a natural experiment is in 
progress in our setting, and our estimates can help us to 
understand the order of magnitude of costs and effect- 
iveness expected from these 'unplanned' changes. 

Regarding the effectiveness measure, different sources 
of uncertainty in the measurement of the outcome- 
adherence relationship should be carefully considered. 
As the data on the outcome onset were drawn from the 
archive of hospitalised patients, our conclusions only 
concern fractures treated on an inpatients basis. This 
issue is particularly important with respect to vertebral 
fractures, the majority of which are not hospitalised in 
most countries. The inclusion of these would have sig- 
nificantly decreased the estimated cost per fracture-free 
year gained. 



Measurement errors due to inaccuracies in the identi- 
fication of therapeutic and diagnostic codes might also 
have occurred. Furthermore, in our study adherence 
with treatment was derived from drug prescriptions, 
which is the most feasible and widely used method to 
estimate it in large populations. 42 With this method, 
however, it must be assumed that the PDC by a prescrip- 
tion corresponds to the proportion of days of drug use, 
which may not be invariably the case. As exposure mis- 
classification is expected to be independent from the 
disease status, the reduction of fracture risk associated 
with high levels of adherence might be larger than that 
quantified by our study 43 and consequently, the 
adherence-enhancing impact should become more 
cost effective. 

As the allocation to adherence levels was not rando- 
mised in our study, the results may be affected by con- 
founding. Notwithstanding this, we attempted to limit 
confounding by adjusting estimates for available demo- 
graphic, therapeutic and clinical factors, and by external 
adjustment for unmeasured confounders, higher adher- 
ence could still be a surrogate for other unmeasured 
characteristics. For example, evidence exists that patients 
with severe osteoporosis appear to be more adherent 
with BP therapy. 44 45 This suggests that the reduction of 
fracture risk associated with high adherence might be 
larger, and enhanced adherence should become more 
cost effective, than that quantified by our study. Other 
factors might exert stronger confounding effect on the 
relation of interest. For instance, patients with cognitive 
impairment, those with family history of hip fracture and 
current smokers, may have poorer adherence and worse 
outcomes. 

Finally, two further weaknesses need to be empha- 
sised. One, in the current application, our study esti- 
mates only concern the time span along which data are 
available (ie, the 6-year period of follow-up). Two, our 
implicit assumption of steady effect of the adherence 
over the entire period of follow-up may be unrealistic. 
This suggests that additional methodological work is 
needed for refining our approach. 

All these limitations taken together, indicate that cost- 
effectiveness profile obtained by our approach may lead 
to overestimation (not capturing outpatient fractures 
and those who are adherent likely having high risk of 
fracture at baseline), as well as to underestimation of 
costs (mainly those of intervention to improve adher- 
ence). For this reason, in the current form our approach 
only indicates the direction, rather than the exact mag- 
nitude, of trends in costs and effectiveness along scen- 
arios of enhanced adherence. 

In conclusion, the present study suggests that our 
approach based on HCU databases is useful to explore 
the cost-effectiveness of enhancing adherence with 
chronic therapies. We found that the impact of 
enhanced adherence offers important benefits in redu- 
cing the fracture risk but at a substantial cost. However, 
proper treatment strategies (eg, increasing adherence 
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monitoring of women at high risk of fracture and dis- 
pensing drugs with cheapest costs), are likely to be asso- 
ciated with added value for the healthcare system. 
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